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Abstract 

A  series  of  optical  depth  calculations  were 
accomplished  to  assess  the  effects  of  various  existing  dust 
and  soot  models  on  the  transmission  of  incident  sunlight. 

A  change  in  the  standard  deviation  of  the  particle  size 
distribution  from  two  to  four,  assuming  constant  total 
density,  resulted  in  a  decrease  in  the  visible  optical  depth 
by  a  factor  of  ten. 

A  technique  using  a  method  of  direct  integration  was 
developed  for  the  calculation  of  the  effective  optical  depth 
of  nuclear  induced  dust  and  soot  clouds.  Contributions  from 
directly  transmitted  photons,  first  scattered  photons  using 
anisotropic  cross-sections,  and  all  subsequently  scattered 
photons  were  used  to  calculate  the  amount  of  light  trans¬ 
mitted  through  the  cloud.  Absorption  effects  were  also 
included. 

The  results  of  this  study  were  comparable  to  the 
results  of  several  recent  nuclear  winter  studies. 


THE  EFFECT  OF  DUST  MODELS  ON 


GLOBAL  NUCLEAR  WINTER 


I.  Introduction 


Background 

Many  scientists  have  expressed  growing  concern  over 
the  climatic  effects  of  a  nuclear  war.  Large  amounts  of 
soot  and  dust  injected  into  the  upper  atmosphere  may  absorb 
and  reflect  large  enough  amounts  of  solar  radiation  for  a 
prolonged  enough  period  of  time  that  the  surface  tempera¬ 
ture  of  the  planet  could  be  affected.  The  low  temperature 
along  with  toxic  smog  and  radioactivity  following  a  nuclear 
war  have  become  known  as  "nuclear  winter"  (Turco  et  al., 
1983a: 33)  . 

Past  studies  in  this  area  (NRC,  1985;  Ramaswamy  and 
Kiehl,  1985;  Turco  et  al.,  1983a)  have  calculated  the  light 
levels  and  temperatures  following  a  nuclear  exchange  for  a 
variety  of  scenarios.  Some  of  the  more  extreme  scenarios 
predict  light  levels  below  the  limits  of  photosynthesis 
and  temperature  drops  of  up  to  30°  C. 

The  optical  properties  of  a  nuclear  cloud  depend  on 
a  variety  of  factors.  For  example,  total  dust  mass,  index 
of  refraction,  composition,  total  area  covered,  height  and 
thickness  of  the  cloud  are  a  few  of  the  important  factors 
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which  need  to  be  determined.  Many  generally  accepted  values 
for  these  parameters  exist  and  have  been  used  in  the  pre¬ 
vious  studies.  But  the  one  factor  which  has  received  little 
or  no  consideration  is  the  sensitivity  of  the  particle  size 
distribution  of  the  cloud  on  the  nuclear  winter  effect. 

Purpose 

The  purpose  of  this  study  is  to  examine  the  effect 
of  particle  size  distribution  and  cloud  composition  on  the 
attenuation  of  solar  radiation  due  to  a  nuclear  cloud. 

Approach 

The  cloud  was  assumed  to  be  a  semi-infinite  cloud 
instantaneously  distributed  over  the  region  of  interest. 
Vertically  the  cloud  is  modeled  as  a  series  of  horizontal 
layers  where  each  individual  layer  is  homogeneous  and  con¬ 
tains  the  appropriate  fraction  of  the  total  mass  lofted. 

The  Air  Force  Geophysics  Laboratory's  Mie  Scattering  com¬ 
puter  code,  MIE3  (Blatner,  1972),  was  used  to  calculate 
the  absorption  and  scattering  coefficients  of  the  nuclear 
cloud.  This  code  was  modified  slightly  to  enable  it  to  run 
on  the  ASD  Cyber  computer  here  at  WPAFB.  Further  modifica¬ 
tions  were  incorporated  to  produce  an  abbreviated  output 
file  containing  a  table  of  scattering  and  extinction  coeffi¬ 
cients  for  the  different  layers  of  the  cloud.  These  values 
were  then  used  by  the  computer  code  SCAT  (see  the  appendix) 
to  calculate  the  fraction  of  visible  sunlight  transmitted 
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through  the  cloud.  This  code  uses  a  method  of  direct  inte¬ 
gration  which  accounts  for  anisotropic  and  multiple  scatter¬ 
ing. 

An  alternative  method  of  calculating  optical  coeffi¬ 
cients  was  also  investigated.  The  Mie  scattering  subroutine 
of  MIE3  was  replaced  with  a  subroutine  using  Diermendjian' s 
formulas  for  extinction  (Diermendjian,  1969:29).  The  accu¬ 
racy  of  this  approximation  was  then  determined  by  comparison 
to  the  results  obtained  using  Mie  theory. 

Overview 

Chapter  II  of  this  report  contains  a  brief  discus¬ 
sion  of  extinction  and  optical  cross-sections.  The  method 
of  direct  integration  is  developed  in  Chapter  III.  Chap¬ 
ter  IV  describes  the  different  types  of  dust  and  soot  clouds 
used  for  the  calculations.  The  results  are  listed  in 
Chapter  V  and  Chapter  VI  is  a  discussion  of  the  results. 

A  summary  and  list  of  recommendations  are  contained  in 
Chapter  VII.  Included  as  an  appendix  is  a  description  and 
listing  of  the  computer  code  SCAT. 
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II.  Theory 


Extinction 

The  interaction  of  sunlight  as  it  passes  through  a 
nuclear  cloud  consists  of  two  processes,  scattering  and 
absorption.  The  two  processes  are  additive  and  their  sum  is 
referred  to  as  extinction. 

Scattering.  When  incident  light  is  scattered  by  a 
dust  particle,  the  particle  becomes  a  point  source  of  the 
scattered  energy.  This  energy  is  reradiated  into  the  total 
solid  angle  around  the  particle.  A  photon  can  therefore  be 
scattered  in  any  direction  and  the  scattering  continues 
until  the  photon  is  lost  out  of  the  top  or  bottom  of  the 
cloud  or  until  it  is  absorbed.  The  light  reaching  the 
ground  is  therefore  a  combination  of  the  directly  trans¬ 
mitted  light,  neither  scattered  nor  absorbed,  and  the  down¬ 
ward  scattered  light.  This  type  of  transmission  is  referred 
to  as  diffuse  transmission  (McCartney,  1976:33). 

Absorption.  The  absorptive  properties  of  a  cloud 
are  determined  by  the  refractive  index  of  the  medium.  The 
index  of  refraction  consists  of  a  real  and  imaginary  part, 
sometimes  referred  to  as  the  scattering  index  and  absorption 
index  respectively  (McCartney,  1976:225).  A  highly  absorp¬ 
tive  particle  will  have  a  large  imaginary  index,  such  as 
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0.66  which  is  typical  for  carbon.  A  less  absorptive  par¬ 
ticle,  such  as  soil,  will  have  a  much  smaller  absorptive 
index  such  as  0.001  (Ackerman  and  Toon,  1981:3362). 


Cross-Sections 

Mie  theory,  developed  by  Gustav  Mie  in  1908,  and 
Diermend j ian 1 s  approximation  provide  two  means  of  computing 
the  optical  cross-section  of  a  particle.  An  efficiency  fac 
tor,  Q,  is  computed  which  relates  the  geometric  cross- 
section  of  a  particle  to  its  optical  cross-section 
(McCartney,  1976:221).  The  microscopic  cross-section  is 
then  found  by  the  relation 

a  =  Qirr2  (1) 


where 

o  =  microscopic  cross-section 
r  =  particle  radius 


For  a  polydispersion  such  as  a  nuclear  cloud  the  microscopic 
cross-sections  are  integrated  across  the  particle  size  dis¬ 
tribution  to  find  the  average  macroscopic  cross-section,  6, 
of  the  cloud.  Using  Equation  (1)  we  obtain 


-  r.Jl 


r  Qn  (r)  dr 


(2) 


where 

8  =  macroscopic  cross-section 
n(r)  =  particle  size  distribution 
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The  computer  code  MIE3  performs  these  calculations 
and  provides  macroscopic  cross-sections  for  extinction, 
scattering  and  absorption.  These  cross-sections  are  more 
commonly  known  as  the  optical  coefficients  8  ,  8„»  and  6 
respectively  and  have  units  of  inverse  length. 

Angular  Coefficients.  The  coefficients  computed  by 
the  above  method  measure  the  ability  of  the  medium  to 
scatter  and  absorb  light  but  they  do  not  indicate  the  direc¬ 
tion  in  which  the  light  will  scatter.  Mie  scattering 
exhibits  several  maxima  and  minima  in  the  scattered  inten¬ 
sity  as  a  function  of  scattering  angle.  The  positions  of 
these  maxima  depend  upon  particle  size  and  wavelength  of 
the  incident  light  (Scott,  1982:1021). 

Angular  scattering  coefficients  may  be  obtained  by 
introducing  the  angular  distribution  function  f(0),  which 
is  the  fraction  of  scattered  radiation  directed  into  a  unit 
solid  angle  in  the  direction  of  6  (Chu  et  al.,  1975: ix)  . 

The  angular  distribution  functions  can  be  represented  in 
the  following  form, 


f  (0) 


ay 

4tt  Li 


a  P  (cos  0) 
n  n 


n=0 


(3) 


where 

PR  =  Legendre  polynomials, 

aR  =  angular  distribution  coefficients,  and 
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6  =  planar  scattering  angle  measured  from  an 
extension  of  the  initial  direction 

The  angular  distribution  coefficients  are  a  function  of  wave¬ 
length,  particle  size  and  index  of  refraction  and  are  listed 
in  Chu  et  al.  (1957).  The  infinite  sum  in  Equation  (3)  can 
usually  be  terminated  after  10  to  15  terms  for  most  of  the 
particle  sizes  of  interest  in  a  nuclear  cloud.  Representa¬ 
tive  values  of  f(e),  from  Equation  (3),  are  listed  in 
Table  I  to  show  the  scattering  pattern  of  visible  light 
(X  =  .5  pm)  from  dust  particles  of  three  different  radii. 

Combining  f(6)  with  8  gives  an  angular  scattering 

s 

coefficient  b  (0)  where 
s 

6s(e)  =  f<e)es  (4) 

Two  limitations  of  this  method  may  be  important 
depending  on  the  composition  of  the  cloud. 

1.  Unless  a  large  data  base  is  used  and  values  of 
f (6)  are  integrated  over  the  entire  size  distribution,  an 
average  f ( 0 )  must  be  computed  from  some  average  particle 
size.  This  may  be  acceptable  for  a  narrow  size  distribu¬ 
tion  but  should  be  used  with  caution  on  a  wider  size  dis¬ 
tribution. 

2.  The  coefficients,  an,  in  Equation  (3)  are  based 
on  non-absorbing  particles.  They  can  be  used  with  good 
accuracy  when  the  extinction  index  is  less  than  0.001  but 
should  not  be  used  with  highly  absorptive  materials  such  as 
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TABLE  I 


ANGULAR  DISTRIBUTION  FUNCTIONS 


i 

Angle 

(Degrees) 

f  (6) 

r=.  1 

r= .  25 

r=l 

i 

0 

.18 

.81 

5.75 

10 

.18 

.75 

2.60 

20 

.17 

.59 

.17 

K 

30 

.15 

.40 

.08 

40 

.13 

.23 

.08 

I 

50 

.11 

.11 

.06 

60 

.09 

.04 

.04 

f 

70 

.08 

.02 

.04 

r 

80 

.07 

.01 

.02 

f 

100 

.05 

.01 

.02 

| 

120 

.05 

.01 

.01 

140 

.06 

.01 

.01 

( 

160 

.07 

.01 

.06 

I 

180 

.07 

.01 

.05 

I 

I 

I 

I 


soot  (Chu  et  al.,  1957:ix).  When  calculations  involving 
soot  are  desired,  using  f  (6 )  =  1/4tt  is  a  good  estimate  since 
soot  particles  have  radii  ranging  from  0.01  to  0.1  pm  and 
symmetric  scattering  can  be  assumed  (McCartney,  1976:21? 
Broyles,  1985:330). 

Optical  Depth 

Optical  depth,  x  ,  is  a  nondimensional  number  used 
to  describe  the  cumulative  depletion  that  a  beam  of  light 
experiences  as  it  travels  through  a  medium  (Pittock  et  al., 
1985:309).  In  terms  of  the  optical  coefficient  of  a 
homogeneous  medium,  the  optical  depth  is 

t  =  8z  (5) 

where  z  is  the  path  length.  The  transmittance  is  then  cal¬ 
culated  using 

T  =  exp (-t )  (6) 

An  optical  depth  of  one  then  reduces  the  beam  to  e”^  (i.e., 

37  percent)  of  its  original  value. 

The  term  optical  depth,  sometimes  called  direct  or 

extinction  optical  depth,  is  used  to  refer  to  removal  by 

both  scatter  and  absorption.  Other  form3  of  optical  depth 

are  often  used  and  refer  to  the  process  by  which  the  light 

is  removed.  These  include  scatter  optical  depth,  t  ,  and 

s 

absorption  optical  depth,  x  .  The  form  which  is  most  often 

cl 

used  in  nuclear  winter  calculations  is  the  effective  optical 
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depth,  igff  This  form  accounts  for  all  the  light  which 
eventually  finds  its  way  out  of  the  bottom  of  the  cloud. 

It  includes  both  directly  transmitted  and  downward  scattered 
light. 


III.  Method 


The  method  of  direct  integration  used  in  this  study 
to  calculate  the  effective  optical  depth  of  a  nuclear  cloud, 
when  scattering  is  significant,  is  adapted  from  the  work 
done  by  Robertson  (1968).  Robertson  demonstrated  that  the 
accuracy  of  this  method  approaches  that  of  Monte  Carlo  cal¬ 
culations  for  the  transport  of  X-ray  photons  in  matter. 

The  procedure  used  by  Robertson  is  outlined  here  along  with 
the  modifications  made  to  adapt  it  to  visible  light  photons 
in  dust  cloud  calculations. 

Geometry 

Figure  1  shows  the  geometry  of  the  semi- infinite 
dust  cloud.  An  arbitrary  target  column  is  selected  and 
divided  into  equally  sized  mesh  cells  of  thickness  Az  and 
diameter  Az.  Each  mesh  cell  of  the  target  column  is 
encircled  by  an  infinite  number  of  concentric  tori  having 

2 

a  radial  thickness  of  Az  and  a  cross-sectional  area  of  Az  . 
Each  torus  acts  as  a  source  kernel  which  will  emit  scattered 
photons  toward  the  target  column.  Field  points  are  located 
at  the  geometric  center  of  each  cell  and  will  serve  as 
reference  points  for  the  calculation  of  photon  fluence.  A 
field  point  is  also  placed  at  the  bottom  of  the  target 
column;  the  fluence  at  this  point  will  be  used  to  calculate 
the  effective  optical  depth. 
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Cloud  Cross-section 


Target  Column 


I 


Assumptions 

1.  The  photon  fluence  at  any  field  point  is  the 


I 

I 

I 

I 

J 

j 

i 

l 

I 

I 

I 

I 


same  as  the  mean  fluence  across  the  cell. 

2.  Anisotropic  scatter  is  used  on  the  first 
scatter  only.  All  subsequent  scatters  are  assumed  isotropic. 

3.  Scatterers  are  spherical  in  shape  and  uniform 
in  density. 

4.  Sunlight  is  assumed  to  have  a  wavelength  of 

0.5  ym. 

5.  Scattered  photons  are  emitted  from  the  point 
sources  located  along  the  toric  center  line  of  each  source 
kernel. 

6.  Extinction  and  scattering  coefficients  remain 
constant  throughout  each  horizontal  subdivision  of  the 
cloud. 

Direct  Fluence 

The  direct  fluence  reaching  a  given  point  along 
the  target  column  is  calculated  using  pure  extinction. 

♦J  =  <P°  exp  [T,  Bj  r.j  (7) 


where 

0  2 
=  direct  fluence  at  point  i  (photons/m  ) 

2 

4>0  =  incident  fluence  (photons/m  ) 
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>•  A. 


ft-1  =  extinction  coefficient  for  mesh  cell  j 
e 

L  =  number  of  mesh  cells  encountered  between  top 
of  cloud  and  point  i 


r 


j 


=  distance  traveled  in  mesh  cell 


j 


The  extinction  optical  depth  can  now  be  calculated 
by  using  the  direct  fluence  reaching  the  last  field  point. 
Using  Equation  (6) 


t  =  ln(T) 


(8) 


where 

T  =  ^bot^O 

<j>bot  =  fluence  at  the  last  field  point 
First  Scatter 

Scattered  photons,  originating  from  the  volume  of 
air  surrounding  the  target  column,  add  to  the  fluence  of  the 
target  column.  The  contribution  from  a  single  source  kernel 
to  the  fluence  at  a  given  field  point  is 


d(*.)  =  .jKBgf  (6) 


(9) 


where 


d(<J>j)  =  fluence  reaching  field  point  j  from  dVol 


<J>^  =  fluence  in  source  kernel  i 
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=  scattering  coefficient  in  source  kernel  i 

f(6)  =  fraction  of  scattered  radiation  directed 
toward  field  point  j 
use  Eq.  (3)  for  first  scatter 
use  1  / 4tt  for  all  other  scatters 

distance  traveled  through  region  k 

extinction  coefficient  for  region  k 

distance  from  source  kernel  i  to  field 
point  j 

dVol  =  volume  of  source  kernel  i 

In  a  semi- infinite  cloud,  the  fluence  in  any  source 
kernel  will  be  identical  to  the  fluence  in  the  correspond¬ 
ing  mesh  cell  of  the  target  column.  Therefore  the  direct 
fluence  calculated  at  each  of  the  field  points  of  the  target 
column  become  the  source  terms  for  Equation  (9). 

The  total  contribution  of  first  scattered  photons 
reaching  a  given  field  point  is  found  by  integrating  Equa¬ 
tion  (9)  . 


<t>j  =  f  d(4>.)  (10) 

Vol 

Point  kernel  integration  is  used  starting  with  the  inner¬ 
most  column  of  source  kernels  and  continuing  outward  until 
some  rejection  criteria  is  met.  The  integration  is  then 
repeated  for  each  of  the  field  points  in  the  target  column. 


15 


1 


Second  and  Subsequent  Scatters 


The  fluence  of  first  scattered  photons  reaching  the 
target  column  now  become  the  source  terms  in  Equation  (9) 
for  the  second  scatter  calculation.  The  second  scattered 
fluence  will  become  the  source  for  the  third  and  so  on  until 
the  contribution  becomes  negligible. 

Effective  Optical  Depth 

To  compute  optical  depth  we  must  know  the  number  of 
photons  per  unit  area  arriving  at  the  last  field  point.  A 
source  kernel  which  is  offset  from  the  target  column  does  not 
"see"  the  same  area  as  a  source  directly  over  the  last  field 
point  (Figure  2) .  For  this  reason  a  correction  factor  of 
cose  must  be  applied  to  Equation  (9)  when  computing  the 
fluence  at  the  last  field  point. 

The  effective  optical  depth  of  the  cloud  can  then 
be  calculated  using  the  relation 


where  ^  is  the  sum  of  the  direct  and  all  the  scattered 
fluence  reaching  the  last  field  point  at  the  bottom  of  the 
cloud. 


IV.  Cloud  Models 


The  dust  cloud  chosen  for  this  study  was  suggested 

by  Turco  et  al.  (1983b: 7)  as  a  likely  product  of  a  massive 

nuclear  exchange.  The  total  yield  expended  consisted  of  a 

combination  of  surface,  near-surface  and  low  air  bursts 

12 

totaling  9400  megatons.  A  total  of  1.5  x  10  g  of  dust 

was  lofted  with  a  height  dependent  injection  profile  shown 

8  2 

in  Figure  3.  A  cloud  area  of  1.1  x  10  km  was  assumed 
which  represents  the  region  of  the  earth's  surface  between 
30°  and  70°  north  latitude.  The  optical  properties  of  the 
dust  cloud  were  based  on  a  refractive  index  of  1.5-.001i  and 

3 

a  density  of  2500  kg/m  . 

The  soot  cloud,  in  Figure  3,  results  from  the  large- 

scale  urban  and  forest  fires  ignited  during  the  attack. 

This  cloud  was  used  as  the  base  case  for  the  NRC  (1985:82) 

12 

study  and  contains  a  total  mass  of  0.18  x  10  g  with  an 

8  2 

area  of  2.5  x  10  km  covering  the  entire  northern  hemis¬ 
phere.  The  soot  has  an  index  of  refraction  of  1.5-.li  and 
a  density  of  1000  kg/m'*. 

Four  representative  particle  size  distributions  were 
selected  for  this  study  and  are  listed  in  Table  II.  The 
three  dust  distributions  were  selected  to  test  the  sensi¬ 
tivity  of  optical  depth  to  the  distribution  parameters.  The 
soot  distribution  was  chosen  to  investigate  the  impact  of  a 
highly  absorptive  medium. 
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TABLE  II 


PARAMETERS  FOR  PARTICLE  SIZE  DISTRIBUTIONS 


Name 

r  (yin) 

0 

a 

Remarks 

TTAPS-Dust 

0.25 

2 

4 

Log  Normal  with 

Power  Law  Tail 

(Turco  et  al.,  1983b:21) 

DELFIC-Dust 

0.204 

4 

- 

Log  Normal 
(Conners,  1985:11) 

NCAR-A-Dust 

1.65 

2.2 

- 

Log  Normal 

(Patterson,  1977:2080) 

NRC-Soot 

0.1 

2 

- 

Log  Normal 
(NRC,  1985:82) 

r  =  mean  radius 
o  =  standard  deviation 
a  =  exponent  of  power  tail 


The  three  dust  particle  size  distributions  were  all 

normalized  to  the  same  total  mass  density.  The  total  mass 
12 

of  1.5  x  10  g  was  distributed  as  shown  in  Figure  3  and  the 
mass  density  computed  for  each  cloud  layer.  These  mass 
densities  were  then  used  in  the  three  dust  distributions 
for  the  optical  depth  calculations. 
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IV 


V.  Results 


Optical  depth  calculations,  using  the  cloud  models 
described  in  the  previous  chapter,  were  made  for  each  of  the 
four  particle  size  distributions.  The  incident  light  was 
assumed  to  be  visible  sunlight  with  a  fixed  wavelength  of 
0.5  pm.  The  mesh  cell  size  chosen  for  the  integration  was 
one  kilometer  for  dust  and  one-half  kilometer  for  soot. 

This  represented  a  distance  of  approximately  one-fifth  and 
one-half  of  the  shortest  mean  free  path  respectively. 

Angular  scattering  coefficients,  from  Equation  (3),  were 
calculated  using  the  mean  radius  of  the  particle  size  dis¬ 
tribution. 

Integration  and  scattering  calculations  were  halted 
when  the  contribution  from  the  next  iteration  represented 
less  than  1  percent  of  the  total  fluence.  This  convergence 
criteria  required  integration  out  to  approximately  35 
kilometers  and  included  from  two  to  eight  scatters  for  the 
dust  clouds.  The  more  absorbent  soot  cloud  only  required 
calculations  out  to  seven  kilometers  with  three  scatters. 

Tables  III  and  IV  list  the  results  of  these  calcu¬ 
lations.  Table  III  compares  the  accuracy  of  Diermendjian' s 
approximation  to  that  of  Mie  theory  for  both  extinction  and 
scattering  optical  depth.  The  results  of  the  multiple 
scattering  calculations  used  to  compute  effective  optical 
depth  are  listed  in  Table  IV. 
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TABLE  III 


DIRECT  OPTICAL  DEPTH 


Cloud 

Model 

Removal 

Process 

Optical  Depth 

%  Error 

Mie 

Diermendjian 

TTAPS 

T- 

2.42 

2.34 

3 

Dust 

6 

T 

S 

2.35 

2.30 

2 

DELFIC 

T 

0.24 

0.23 

4 

Dust 

e 

T 

S 

0.21 

0.21 

0 

NCAR-A 

T 

1.00 

0.96 

4 

Dust 

e 

Ts 

0.89 

0.88 

1 

NRC 

T- 

3.94 

3.82 

3 

Soot 

6 

T 

2.52 

2.79 

11 

S 

The  extinction  optical  depth  is  listed  along  with  the  opti¬ 
cal  depth  which  results  after  each  successive  scatter  until 
the  rejection  criteria  is  met.  Table  IV  was  calculated 
using  the  Mie  coefficients.  Diermendjian* s  approximation 
produced  similar  results  with  errors  ranging  from  5  to  7 
percent. 
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TABLE  IV 

EFFECTIVE  OPTICAL  DEPTH 


Number  of 
Scatters 

Particle  Size 

Distribution 

TTAPS 

Dust 

NCAR-A 

Dust 

NRC 

Soot 

DELFIC 

Dust 

Direct 

2.42 

1.00 

3.94 

0.24 

1 

1.70 

0.86 

3.80 

0.12 

2 

1.50 

0.80 

3.73 

0.11 

3 

1.38 

0.77 

3.70 

II 

4 

1.31 

0.76 

II 

II 

5 

1.26 

0.75 

II 

II 

6 

1.23 

II 

•1 

II 

7 

1.21 

II 

II 

II 

8 

1.20 

It 

II 

II 
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VI .  Discussion 


This  chapter  deals  with  four  main  topics:  (1)  the 
characteristics  of  the  cloud  and  how  they  effect  optical 
depth,  (2)  the  validity  of  the  basic  assumptions  used  in 
the  calculations,  (3)  comparison  of  the  results  to  previous 
work,  and  (4)  the  accuracy  of  Diermendjian' s  formulas. 

Cloud  Characteristics 

The  main  characteristic  of  a  dust  cloud  which  is 
altered  by  varying  the  particle  size  distribution  param¬ 
eters  is  the  number  concentration  of  submicron  particles 
(i.e.  0.1  <  r  <  1  ym) .  Submicron  particles  are  extremely 
efficient  scatterers  with  optical  cross-sections  that  can 
be  as  much  as  4.4  times  larger  than  their  geometric  cross- 
sections  (McCartney,  1976:248). 

The  TTAPS  dust  and  NRC  soot  clouds  are  character¬ 
ized  by  very  narrow  size  distributions  with  high  concentra¬ 
tions  of  submicron  particles.  In  comparison  the  DELFIC  dust 
and  NCAR  dust  clouds  use  distributions  with  a  much  lower 
fraction  of  small  particles.  This  lower  concentration 
results  from  the  larger  standard  deviation  of  the  DELFIC 
distribution  and  the  larger  mean  radius  of  the  NCAR  dis¬ 
tribution. 

The  impact  of  these  submicron  particles  can  be  seen 
by  comparing  the  optical  depths  listed  in  Table  III.  The 
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three  dust  clouds  show  the  sensitivity  of  optical  depth  to 
the  distribution  parameters.  Comparison  of  the  TTAPS  and 
DELFIC  dust  results  indicates  that  the  standard  deviation 
is  the  most  sensitive  parameter.  An  increase  in  standard 
deviation  by  a  factor  of  two  decreases  the  optical  depth  by 
a  factor  of  ten.  The  NCAR  dust  distribution  has  approxi¬ 
mately  the  same  standard  deviation  as  the  TTAPS  but  the 
mean  radius  is  larger  by  a  factor  of  six?  this  also 
decreases  the  optical  depth  but  not  to  such  a  great  extent. 
Although  a  direct  comparison  of  the  dust  and  soot  is  not 
possible,  due  to  the  different  refractive  indices,  the 
trend  toward  high  scattering  efficiencies  can  be  seen  in 
both  the  TTAPS  dust  and  NRC  soot  distributions. 

The  results  listed  in  Table  IV  indicate  that  the 
number  of  scatters  which  must  be  considered  is  also 
affected  by  the  cloud  properties.  The  cloud's  scattering 
efficiency  and  absorptive  index  will  determine  how  much 
light  is  scattered  and  the  rate  at  which  this  scattered 
light  is  attenuated.  The  source  terms,  used  to  calculate 
the  scattered  contribution  to  the  total  fluence,  originate 
from  the  fraction  of  light  scattered  out  of  the  direct  beam. 
Both  the  NRC  soot  and  TTAPS  dust  clouds  scatter  approxi¬ 
mately  91  percent  of  the  incident  light  and  therefore  have 
relatively  large  source  terms  compared  to  the  DELFIC  dust 
cloud  which  only  scatters  21  percent  of  the  incident  light. 
The  scattered  light  in  the  soot  cloud  is  quickly  attenuated, 
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after  only  three  scatters,  due  to  the  highly  absorptive 
nature  of  the  soot.  In  the  dust  cloud,  however,  less  than 
1  percent  of  the  light  is  absorbed  and  the  light  continues 
to  scatter  until  spherical  divergence  becomes  a  factor. 

As  a  result  of  this  scattering  the  optical  depth  of  the 
dust  cloud  is  reduced  by  as  much  as  50  percent,  while  in  the 
soot  cloud  the  scattered  light  is  unable  to  make  a  signifi¬ 
cant  contribution  and  only  reduces  the  optical  depth  by 
6  percent. 

Validity  of  Assumptions 

The  basic  assumption  of  a  semi-infinite  cloud 
proved  to  be  very  good.  The  integration  in  the  radial 
direction  converged  in  from  one-half  to  seven  mean  free 
paths.  The  soot  cloud  required  integration  to  seven 
kilometers  which  represents  approximately  seven  mean  free 
paths.  Convergence  in  the  less  turbid  clouds  occurred  in 
less  than  one  mean  free  path  due  to  the  small  amount  of 
scattered  light  combined  with  the  effects  of  spherical 
divergence . 

The  photon  fluence  as  calculated  in  Equation  (7)  was 
a  good  approximation  of  the  mean  fluence  across  the  cell. 

In  all  but  the  soot  cloud  a  cell  thickness  of  one  kilometer 
could  be  used.  The  shorter  mean  free  path  in  the  NRC  soot 
cloud  required  a  cell  thickness  of  one-half  kilometer  to 
keep  the  error  below  1  percent. 
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To  test  the  validity  of  calculating  the  angular 
scattering  coefficients  using  the  mean  radius  of  the  dis¬ 
tribution  an  attempt  was  made  to  bracket  the  correct 
answer.  The  two  extremes  would  be  the  use  of  a  very  small 
radius  and  a  very  large  radius  as  the  basis  for  the  angular 
scattering.  A  small  radius  would  reduce  the  amount  of  for¬ 
ward  scattering  and  tend  to  increase  the  optical  depth 
since  more  light  is  being  scattered  up  and  away  from  the 
source  column.  The  larger  radius  would  also  tend  to 
increase  the  optical  depth,  since  the  extreme  forward  scat¬ 
tering  permits  only  those  source  kernels  directly  adjacent 
to  the  target  column  to  scatter  light  into  the  column. 

Since  both  extremes  yield  optical  depths  greater  than  the 
true  optical  depth  we  can  use  this  as  an  upper  limit.  These 
calculations  were  performed  using  the  TTAPS  dust  cloud  with 
radii  of  0.1  ym  and  1  ym  giving  optical  depths  of  1.39  and 
1.59  respectively.  To  find  the  lower  limit  several  calcula¬ 
tions  were  performed  using  radii  between  0.1  ym  and  1  ym  and 
a  minimum  value  of  1.19  was  found  using  a  radius  of  0.18  ym. 
The  actual  optical  depth  should  therefore  have  a  value 
between  1.19  and  1.39.  Using  the  mean  radius  of  the  size 
distribution  a  value  of  1.20  was  obtained.  Similar  calcula¬ 
tions  were  done  using  the  DELFIC  dust  cloud  and  the  minimum, 
mean  radius  and  maximum  optical  depths  were  0.11,  0.11  and 
0.15  respectively.  The  actual  optical  depth  can  therefore 
be  approximated  with  fairly  good  accuracy  by  using  the  mean 
radius  for  angular  scattering. 
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One  last  assumption  to  be  considered  is  the  use  of 
isotropic  scatter  on  all  scatters  after  the  first.  In  all 
but  the  TTAPS  dust  cloud  the  direct  and  first  scattered 
fluence  accounted  for  over  90  percent  of  the  transmitted 
fluence.  Considering  the  small  contribution  from  second 
and  subsequent  scatters  and  the  difficulty  in  tracking  the 
direction  of  each  scatter,  the  use  of  isotropic  scattering 
is  not  a  bad  assumption. 

Previous  Work 

Both  the  TTAPS  dust  and  NRC  soot  clouds  were  used 
as  the  basis  for  nuclear  winter  calculations  in  previous 
studies.  Turco  et  al.  (1983b)  used  the  TTAPS  cloud  as  a 
baseline  scenario  and  reported  an  effective  optical  depth 
of  1.4,  whereas  1.20  was  found  here.  A  radiative  trans¬ 
fer  model,  developed  by  Sagan  and  Pollock  (1976),  was  used 
as  an  approximation  to  account  for  the  combined  effects  of 
multiple  scattering  and  absorption.  Their  model  does  not 
use  specific  scattering  angles  but  instead  uses  a  single 
factor  to  account  for  the  fraction  of  light  scattered  in  the 
forward  direction.  The  different  approach  to  angular 
scattering  used  in  this  study  may  account  for  the  slight 
variation  in  results. 

The  NRC  study  reported  an  extinction  optical  depth 
of  4.0  for  their  baseline  smoke  cloud  (NRC,  1985:84).  This 
compares  favorably  with  the  extinction  optical  depth  of 
3.94  that  was  calculated  in  this  study. 


Diermendjian'  s  Approximation 

Diermendjian' s  formulas  provide  a  simple  alternative 
means  of  calculating  optical  efficiencies.  The  results 
listed  in  Table  III  show  that  this  method  is  accurate  to 
within  4  percent  for  most  cases.  Larger  errors  have  been 
reported  (Broyles,  1985:328)  when  the  index  of  refraction 
exceeds  1.5-.25i.  This  may  become  a  factor  for  soot  and 
smoke  cloud  calculations  which  typically  involve  imaginary 
indices  of  0.1  to  0.6. 


II.  Summary  and  Recommendations 


A  method  of  computing  the  effective  optical  depth 
of  a  nuclear  cloud  was  developed.  The  method  uses  direct 
integration,  anisotropic  and  multiple  scattering  to  deter¬ 
mine  the  fraction  of  incident  light  which  emerges  from  the 


bottom  of  the  cloud.  The  results  of  this  method  are 


comparable  to  the  results  obtained  in  previous  nuclear 
winter  studies. 


A  comparative  study  was  performed  to  determine  the 
effect  of  particle  size  distribution  parameters  on  the 
nuclear  winter  phenomenon.  It  was  found  that  the  standard 
deviation  of  the  distribution  was  the  most  sensitive  factor. 


An  increase  in  the  standard  deviation  by  a  factor  of  two 
can  decrease  the  optical  depth  by  a  factor  of  ten. 

The  relationship  between  the  impact  of  multiple 
scattering  and  cloud  composition  was  also  investigated. 

The  absorptive  properties  of  the  cloud  played  an  important 
role  in  determining  the  reduction  in  optical  depth  due  to 
multiple  scattering.  For  a  nonabsorptive  dust  cloud  the 
optical  depth  can  be  reduced  by  as  much  as  50  percent,  but 
in  a  more  absorptive  soot  cloud  the  reduction  was  only  6 
percent. 


Recommendations 


Future  investigations  of  the  nuclear  winter 
phenomenon  should  include  a  study  of  the  time  dependent 
aspects  of  optical  depth.  These  should  include  sedimenta¬ 
tion,  agglomeration,  rainout  and  horizontal  diffusion. 
Diermendjian' s  formulas  could  be  used  in  place  of  the 
lengthy  Mie  scattering  code,  and  provide  a  simple  and 
fairly  accurate  means  of  recalculating  optical  coefficients 
as  the  cloud  properties  change  with  time. 

To  reduce  the  overall  uncertainties  in  nuclear 
winter  calculations,  the  physical  and  optical  properties 
of  dust  and  soot  should  be  more  accurately  defined. 
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J  Appendix:  Computer  Program  Listing  (SCAT) 

|  The  following  pages  contain  a  listing  of  the  com¬ 

puter  program  SCAT  used  for  the  calculation  of  the  effec- 

|  tive  optical  depth  of  a  nuclear  cloud.  The  program  is 

written  in  FORTRAN  IV  and  will  run  on  the  ASD  CYBER  com¬ 
puter. 
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PROGRAM  SCAT  (INPUT, OUTPUT, TAPE9, TAPE7 ,TAPE5*INPUT,TAPE6*OUTPUT) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


DIMENSION  BEX  (20) , BSC  (20) , BE (100) ,BS(100) ,F0(100) , FI (100) 

PROGRAM  COMPUTES  EFFECTIVE  OPTICAL  DEPTH  OF  AN  AEROSOL  CLOUD 
USES  ANISOTROPIC  SCATTER  ON  FIRST  SCATTER  ONLY  ALL  OTHER 
SCATTERS  ARE  ASSUMED  ISOTROPIC 

USER  SETS  LIMITS  OF  INTEGRATION,  INTEGRATION  INCREMENT  AND 
NUMBER  OF  SCATTERS 

VARIABLES  DEFINED  BY  USER 

KTHICK :  CLOUD  THICKNESS  (KM) 

KSLAB :  THICKNESS  OF  CLOUD  LAYERS  (KM) 

NDIST:  RADIAL  DISTANCE  TO  WHICH  INTEGRATION  IS  TO 
BE  ACCOMPLISHED  (KM)  ,DISTANCE=NDIST*KDZ 
NSCAT:  NUMBER  OF  SCATTERS  TO  BE  USED 
KDZ :  INTEGRATION  INCREMENT  (KM) 

INPUT  FROM  TAPE9 

TAPE 9  SHOULD  CONTAIN  ( KTH ICK/KSLAB )  SETS  OF 
EXTINCTION  AND  SCATTERING  COEFFICIENTS  (KM**-1) 
F0RMAT(2E15.5) 

DERIVED  VARIABLES 

NCOEFF :  NUMBER  OF  COEFFICIENT  SETS  TO  BE  READ  FROM  TAPE9 
NDZ:  NUMBER  OF  INTERNAL  FIELD  POINTS 
KBOT:  NUMBER  TO  DEFINE  LAST  FIELD  POINT 
COSIN:  COSINE  OF  SCATTERING  ANGLE 

R:  DISTANCE  BETWEEN  SOURCE  AND  FIELD  POINT 
FTOT:  TOTAL  FLUENCE  AT  LAST  FIELD  POINT 

OUTPUT  IS  WRITTEN  TO  TAPE 7 

**  DEFINE  CLOUD  ** 

KTHCK=28 

KSLAB=2 

KDZ=1 

NDIST=35 

NSCAT=11 

NCOEFF =KTH ICK/KSLAB 
NDZ=KTH ICK/KDZ 
KBOT=NDZ+l 

**  READ  COEFFICIENTS  FROM  TAPE 9  ** 

READ(9, 100)  (BEX (J) ,BSC(J) , J-l, NCOEFF) 


33 


oo  n  on  o  o  oo 


**  DIVIDE  CLOUD  INTO  NDZ  LAYERS  AND  ASSIGN  EXTINCTION** 
**  AND  SCATTERING  COEFFICIENTS  TO  EACH  LAYER  ** 

INC*KSLAB/KDZ 
1=0 
K=INC 

DO  20  L=1 , NCOEFF 
II=I+L 
KK=K+L 

DO  10  N=II , KK 

BE  (N)  =BEX  (L) 

BS (N) =BSC (L) 

CONTINUE 
1=1+ INC 
K=K+INC 
CONTINUE 

**  INITIALIZE  ARRAYS  ** 

DO  30  N=1 ,KBOT 
F0 (N) =0.0 
FI (N) =0.0 
CONTINUE 
BE(KBOT) =0.0 
BS (KBOT) =0 . 0 

**  COMPUTE  DIRECT  FLUENCE  ** 

F0 (1) =EXP (-BE (1) *KDZ/2 . ) 

DO  40  N=2 ,NDZ 

F0  (N)  =F0  (N-l)  *EXP("(BE  (N-l)  +BE  (N) )  *KDZ/2 . ) 

CONTINUE 

FO (KBOT) =F0 (KBOT-1) *EXP  (-BE (KBOT-1) *KDZ/2 . ) 

FTOT=FO (KBOT) 

WRITE (7, 200)  NDIST*KDZ,KDZ 
WRITE (7, 300) 

WRITE (7, 400) 

WRITE (7 , 500)  FO  (KBOT)  ,ALOG(FO  (KBOT) ) 

**  INTEGRATE  OVER  VOLUME  OF  CLOUD  ** 

DO  85  NS=1,NSCAT 
DO  75  L=1 ,NDIST 

DO  70  N=l,KBOT 

IF  (N  .EQ.  KBOT)  A=KDZ/2. 

IF  (N  .NE.  KBOT)  A=0.0 
DO  65  M=1,NDZ 

R=SQRT( (L*KDZ) **2+(  (M-N-A)  *KDZ)  **2) 

COSIN- ( (N-M-A) *KDZ)  /R 
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**  COMPUTE  EXTINCTION  (R2)  BETWEEN  SOURCE  M  AND  CELL  N  ** 
IF  (M  .EQ.  N)  R2«BE(N)*KDZ*L 
IF  (M  .EQ.  N)  GO  TO  60 
R2-0.0 

IF  (M  .GT.  N)  GO  TO  50 
DO  45  K“M,N 

R1 “FLOAT (KDZ) 

IF  (K  .EQ.  M  .OR.  K  .EQ.  N)  Rl=KDZ/2. 

R2=R2+R1*BE (K) /ABS (COSIN) 

45  CONTINUE 

GO  TO  60 

50  DO  55  K=N,M 

R 1 “FLOAT (KDZ) 

IF  (K  .EQ.  N  .OR.  K  .EQ.  M)  Rl=KDZ/2. 

R2=R2+R1*BE (K) /ABS (COSIN) 

55  CONTINUE 

60  CONTINUE 

**  COMPUTE  SCATTERED  FLUENCE  ARRIVING  AT  CELL  N  ** 

F*F0  (M)  *BS (M)  *L*KDZ**3*EXP  (-R2)  /  (2 .  *R**2) 

IF  (N  .EQ.  KBOT)  F=F*COSIN 
IF  (NS  .EQ.  1)  F=F*FRAC (COSIN) 

F1(N)=F1(N)+F 
65  CONTINUE 

70  CONTINUE 
75  CONTINUE 

FTOT=FTOT+F  1  (KBOT) 

WRITE (7, 600)  Fl (KBOT) ,NS,AL0G (FTOT) 

DO  80  N=1 , KBOT 
FO  (N)  =F1  (N) 

Fl (N)“0.0 
80  CONTINUE 
85  CONTINUE 

WRITE (7 ,400) 

WRITE (7, 700)  FTOT 
100  FORMAT (2E15. 5) 

200  FORMAT (15 , *  KM*,8X,*DZ=*,I3/) 

300  FORMAT  (10X,*%FLUX*,4X,*  SOURCE*,  8X,*  OPTICAL  DEPTH*) 

400  FORMAT  (9X,* - *) 

500  FORMAT  (F15 . 5 , 5X ,  ‘DIRECT* ,  3X  ,F15 . 5) 

600  FORMAT  (F15 . 5 , 2X ,  13  ,*  SCATTER* ,  F15 . 5) 

700  FORMAT (F15. 5 ,4X, ‘TOTAL*) 

STOP 

END 
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FUNCTION  FRAC (COSIN) 

DIMENSION  A(26) ,P(26) 

C 

C  COMPUTES  ANGULAR  DISTRIBUTION  FUNCTION 

C  COSIN:  COSINE  OF  SCATTERING  ANGLE 

C  P:  LEGENDRE  POLYNOMIALS 

C  A:  ANGULAR  DISTRIBUTION  COEFFICIENTS 

C  NA:  THE  NUMBER  OF  COEFFICIENTS  BEING  USED 

C  ALPHA:  2*PI*R/WAVELENGTH 

C 

C  **  DATA  FOR  ALPHA*  1  USE  NA=5  ** 

C  DATA  A/1.0,. 59685,. 56415,. 10955,. 00902/ 

C 

C  **  DATA  FOR  ALPHA* 3  USE  NA=11  ** 

C  DATA  A/1.0,2.203,2.58189,2.18461,1.34801,.56111,.21739, 

C  $  .04954,. 0083,. 00107,. 00012/ 

C 

C  **  DATA  FOR  ALPHA*10  USE  NA=26  ** 

DATA  A/1.0,2.22873,3.38409,3.78255,4.76483,5.16916,5.70666, 
$  5.89478,5.92665,5.875,5.41723,5.11591,4.30744, 

$  3. 77023, 2. 78576, 2. 14168, 1.13952,. 86443,. 61407, 

$  . 61738 ,. 76154 , . 48723 , . 62348 , - . 09834 , . 00737 , . 00203/ 

C 

NA=26 
FRAC=0 . 0 
P(l)-1.0 
P (2) =COSIN 
DO  20  N=3,NA 

P (N) *(P(N-1) * (2*N-3) *COSIN-P (N-2) * (N-2) ) / (N-l) 

20  CONTINUE 

DO  30  N=1 ,NA 

FRAC=FRAC+A(N)*P(N) 

30  CONTINUE 
RETURN 
END 
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